The transient and the late time attractor tachyon dark energy: Can we distinguish it 

from quintessence ? 
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The string inspired tachyon field can serve as a candidate of dark energy. Its equation of state 
parameter w varies from to —1. In case of tachyon field potential V(<j>) — > slower(faster) than 
l/4> 2 at infinity, dark energy (dark matter) is a late time attractor. We investigate the tachyon dark 
energy models under the assumption that w is close to —1. We find that all the models exhibit 
unique behavior around the present epoch which is exactly same as that of the thawing quintessence. 

PACS numbers: 



I. INTRODUCTION 

One of the most challenging problems of modern cos- 
mology is associated with late time acceleration of uni- 
verse which is supported by observations of complemen- 
tary nature. According to the standard lore, an ex- 
otic perfect barotropic fluid with large negative pres- 
sure dubbed dark energy can account for repulsive ef- 
fect causing acceleration[l|, 0, H[ ■ The simplest example 
of dark energy is provided by cosmological constant A. 
The model is consistent with observation but is plagued 
with difficult theoretical issues. The field theoretic un- 
derstanding of A is far from being satisfactory and its 
small numerical values gives rise to problems of fine tun- 
ing and coincidence. A variety of scalar field models in- 
cluding quintessence, tachyons, phantoms and K-csscncc 
has been investigated in the recent years to address the 
problem @, i, These models have some advantage 
over the cosmological constant: (i) They can mimic cos- 
mological constant at the present epoch and can give rise 
to other observed values of the equation of state param- 
eter w (recent data indicate that w lies in a narrow strip 
around w = w\ = — 1 and is consistent with being below 
this value), (ii) They can alleviate the fine tuning and 
coincidence problems. 

The scalar field model, which is the simplest gener- 
alization of cosmological constant, is one with a linear 
potential @ . This model starts with a cosmological con- 
stant like behaviour where the scalar field is frozen ini- 
tially due to Hubble damping. Later on, it starts rolling, 
but because the potential has no minimum, it leads to 
a collapsing universe in future. Hence universe in this 
model, has a finite history. 

The more complicated scalar field models can broadly 
be classified into two categories. Models in which scalar 
field mimics the background (radiation/matter) being 
subdominant for most of the evolution history. Only 
at late times it becomes dominant and accounts for the 
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late time acceleration. Such a solution is referred to as 
tracker. In this case w{<j)) ~ Wb = 0, 1/3) before the 
transition from matter like regime or scaling regime to 
accelerated expansion. Tracker models are independent 
of initial conditions used for field evolution but do re- 
quire the tuning of the slope of the scalar field potential. 
During the scaling regime the field energy density is of 
the same order of magnitude as the background energy 
density. 

In second class of models, trackers are absent. Hence 
at early times, the field gets locked (w((f>) = —1) due to 
large Hubble damping and waits for the matter energy 
density to become comparable to field energy density 
which is made to happen at late times. The field then be- 
gins to evolve towards larger values of w((j)) starting from 
w(<j)) = — 1. In this case, for a viable cosmic evolution, 
one chooses ~ p\ during the locking regime which 
requires the tuning of initial conditions of the field. The 
two classes of scalar fields are called Freezing and Thaw- 
ing models. 

In case of standard scalar field (quintessence), there is 
a variety of models which possess tracker solutions. In 
case of tachyon field @, @] (motivated by string theory), 
there exists no solution which can mimic scaling mat- 
ter/radiation regime [{J M, 03, OJ, 0, HI El, 03 E3- 
These models necessarily belong to the class of thawing 
models. Tachyon models do admit scaling solution in 
presence of a hypothetical barotropic fluid with negative 
equation of state. Tachyon fields can be classified by the 
asymptotic behavior of their potentials for large values 
of the field: (i) V((f)) — > faster than l/(/> 2 for — > oo. 
In this case dark matter like solution is a late time at- 
tractor. Dark energy may arise in this case as a transient 
phenomenon, (ii) V(4>) — > slower than 1/0 2 for (j> — > oo 
; these models give rise to dark energy as late time at- 
tractor. The two classes are separated by V(<p) ~ 1/0 2 
which is scaling potential with w(<f>) — const. 

Since observationally, the equation of state parame- 
ter of dark energy is very close to one, we can use this 
information to simplify the dynamics. In case of thaw- 
ing quintessence and phantom field, it allows to obtain 
a generic expression for w which represents the entire 
class of quintessence and phantom models [HI, [l9j]. In 
this paper we apply the same technique to tachyon field 



2 



which belongs to the class of thawing models. With the 
current state of observation, we address the issue of dis- 
tinguishing the tachyon dark energy from the case of 
quintessence. 



II. DYNAMICS OF TACHYON FIELD 

In what follows we shall be interested in the cosmolog- 
ical dynamics of tachyon field which is specified by the 
Dirac-Born-Infcld (DBI) type of action given by 



(1) 



where on phenomcnological grounds, we shall consider 
a wider class of potentials satisfying the restriction that 
V((j)) — > as — > oo. The parameter e = ±1 where 
the plus sign corresponds to the normal tachyon field 
which is non-phantom whereas with minus sign, one can 
model phantom type tachyon fields phcnomenologically. 
In FRW background, the pressure and energy density of 
<f> are given by 



P* = -V(<l>)<Jl-e<tr 



V(<j>) 



P<t> 



The equation of motion which follows from ([I]) is 
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where H is the Hubble parameter 
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The evolution equation can be cast in the following au- 
tonomous form for the convenient use 
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FIG. 1: Plot of dark energy equation of state parameter w$ 
versus £1$ for < f2^ < 0.7 in case of different values of Ao 
for nonphantom case i.e e = 1. The curves are for the po- 
tentials V(4>) — (f)" 3 (dotdashed curve), V(<f>) = <f>~ 2 (dashed 
curve), V((f>) = <ji~ 1 (dotted curve). The black solid line is for 
our analytical approximation (17). The sets (a), (b) and (c) 
are for Ao = 1, 2/3, 1/2 respectively. 



An important remark on the autonomous system is in 
order. Let us consider the inverse power law type poten- 
tial V((j)) = V (j) n (n < 0). Eq.© tells us that T > 3/2 
if n < —2 allowing A to increase monotonously for large 
values of the field. In this case <p — > 1 or w — > where as 
w approaches the de-Sitter limit for n > —2 (r < 3/2). 
These two classes of tachyon potentials are separated by 
the inverse square potential with constant A (r = 3/2) 
which provides the analog of scaling potential in case of 
tachyon. However, there is major difference that in the 
present case, the field can only mimic a hypothetical fluid 
with negative equation of state leading to accelerated ex- 
pansion. Unfortunately, the mass scale in the potential 
turns out be larger than the Planck mass. The class of 
potentials designated by — 2 < n < is free from this 
problem and gives rise to dark energy as late time at- 
tractor of dynamics. In the analysis to follow, it will be 
convenient to use the following quantities 
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where prime denotes the derivative with respect to ln(a). 
Here 76 is defined as pb = (76 — l)p6 for the background 
field. In our subsequent calculations, we shall assume a 
non-relativistic matter for our background field for which 

7b = 1- 



where w = — is the equation of state for the tachyon 
field. One can now express the autonomous equations 
through them: 



7 ; = -6741 -e- /4> ) + 2^3^n;\(l-e^)i (11) 
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FIG. 2: Plot of dark energy equation of state parameter w$ 
versus a in case of different values of Ao for non-phantom 
case i.e e = 1. The curves are for the potentials V((j>) — 
<j>~ 3 (dotdashed curve), V(cj>) = (f>~ 2 (dashed curve), V(<p) = 
<j>~ (Dotted curve). The black solid line is for our analytical 
approximation (17) together with (20). The sets (a), (b) and 
(c) are for Ao = 1,2/3, 1/2 respectively. 



A' = - ev /3^A 2 (l-e 7 ^(r--) (13) 

The first two equations can be combined into one by a 
change of variable from a — > f2^ 
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dj<p = 70 = -2 70 (1 - €7^) 2^370» A(1 - ej<f,) * 

n^(i - n^)(i - e 70 ) 3fi^(i - n^)(i - e 70 ) 

(14) 



A. Late time evolution 

From Eq. (fT4")) , one can see that for non-phantom and 
phantom cases, i.e e = ±1, the equation is completely dif- 
ferent and hence one expects to have different evolutions 
for 70(^0) for non-phantom and phantom cases. 

But we are interested in the investigations of cosmolog- 
ical dynamics around the present epoch where 7^ << 1. 
Secondly, in our case w(<j)) improves slightly beginning 
from the locking regime, thereby, telling us that the slope 
of the potential does not change appreciably. This im- 
plies that the potential is very flat around the present 
epoch such that 



FIG. 3: Plot of dark energy equation of state parameter 
versus Q.^ for < Q$ < 0.7 in case of different values of 
Ao for phantom case i.e e — — 1. The curves are for the po- 
tentials V(4>) — </>~ 3 (dotdashed curve), V(<j>) = (j>~ 2 (dashed 
curve), V{cj>) = (j)- 1 (dotted curve). The black solid line is for 
our analytical approximation (17). The sets (a), (b) and (c) 
are for Ao = 1, 2/3, 1/2 respectively. 



Eq. lfTS")) define the slow roll parameters which allows to 
neglect the <f> term in equation of motion for <f>. In the 
present context, unlike the case of inflation, the evolution 
of field begins in the matter dominated regime and even 
today, the contribution of matter is not negligible. The 
traditional slow roll parameters can not be connected to 
the conditions on slope and curvature of potential which 
essentially requires that Hubble expansion is determined 
by the field energy density alone. Thus the slow roll pa- 
rameters are not that useful in case of late time accelera- 
tion, though, Eq. (fT5|) can still be helpful. In view of the 
aforesaid, we can drop all the terms of order higher than 
70 in Eq. (fT4]) and assume that the slope of the potential 
is constant, A = Ao- These follow from the two slow-roll 
conditions (15) as we shall show later. Evolution equa- 
tion then simplifies to 
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In case of field domination regime, the two conditions in 



Let us note that Eq. (|16p is same as its counter part in case 
of quintessence though the full Eq.(Q3]) is different. The 
difference between tachyon and quintessence dynamics 
is represented by terms of higher order than 70. Thus 
if we restrict our investigation of dark energy dynamics 
very close to cosmological constant behavior, we can not 
distinguish tachyon dark energy from quintessence. Also 
equation (16) is independent of e. Hence (1 + w) for 
non-phantom case and —(1 + w) for phantom case, have 
exactly similar evolution around cosmological constant. 
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FIG. 4: Plot of dark energy equation of state parameter 
Wtf, versus a in case of different values of Ao for phantom 
case i.e e = —1. The curves are for the potentials V((j>) = 
<f>~ 3 (dotdashed curve), V(4>) = (\T 2 (dashed curve), V((p) = 
<j>~ (Dotted curve). The black solid line is for our analytical 
approximation (17) together with (20). The sets (a), (b) and 
(c) are for A = 1,2/3, 1/2 respectively. 



poses. One can also arrives the same behavior for phan- 
tom case. In figure 1 and figure 3, we show the our analyt- 
ical approximation for w(ft$) in comparison with the nu- 
merical solutions of the exact equations for different po- 
tentials with different initial values for A for non-phantom 
and phantom cases. They show that our approximation 
works reasonably well as long as Ao is small, i.e as long 
as the slow-roll conditions are satisfied. 

Next, we can use eqn (12) to solve for fl ! j ) (a) to deter- 
mine w(a). assuming 7^ << 1, this gives 



a, 
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00 
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(19) 



where fi^o is the present day value of £7^. Equation (17) 
and (19) gives the complete behavior for the equation 
of state w(a) for tachyon fields with potentials satisfying 
the slow-roll conditions (15). One can also express the 
parameter Ao in terms of the present day value wq of the 
equation of state which is quite straightforward. This 
behaviours are shown in figure 2 and figure 4 for non- 
phantom and phantom cases. 

Similar to the case of thawing quintessence, non- 
phantom tachyon models are restricted to a part of the 
w' — w plane. To specify the the limits, let us define a 
parameter X 



Eq. (|16p can be transformed into a linear differential 
equation with the change of variable s 2 = 7</>, we have 
boundary condition 7^ = at f2^ = 0. The resulting 
solution expressed in terms of w((j>) 
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(17) 



Under the approximation 7^ << 1 which is justified 
about the present epoch, all the tachyon models follow a 
general track irrespective of the particular field potential. 
One can see from (17) that 1+w ~ 0(\ 2 ). Hence the first 
slow roll condition (A << 1) ensures that 1 + w << 1. 
We can quantify our second assumption that the slope 
of the potential does not change appreciably during the 
evolution as A /A << 1. Noting that 7 ~ A 2 and also 
7 << 1, one can then use eqn (13) to write 



V 2 



2 V 3 
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(18) 



together with the first slow-roll condition, this ensures 
the second slow-roll condition to be satisfied. We also 
show in figure 5, the actual behavior of A for different 
potentials for non-phantom case. This also shows A is 
constant during the entire evolution for all practical pur- 



H(pw 



2w(l + w) 



w' = -2Xw(l 



Since the Hubble parameter is determined by matter 
dominated regime in the beginning of evolution, we find 
that X = — 3/2w < 3/2 as w > — 1 which leads to 
the upper limit, w' < 3(1 + w). The lower bound on 
w' is estimated numerically (demanding that at present 
Sl0 <= 0.8) as, if' > —.8(1 + w) giving rise to the per- 
missible region of w'-w plane 



0.8(1 + w) < w' < 3(1 + 10). 



(20) 



In figure 3 we have shown this permissible region to- 
gether with the actual behavior for different potentials. 



III. OBSERVATIONAL CONSTRAINT 

The solution given by Eqs.fTT]) & (jTSJ) for the equa- 
tion of state parameter w versus the scale factor a for 
tachyon field under slow-roll conditions is exactly simi- 
lar to that for a canonical scalar field as obtained ear- 
lier in [l8|, [l||. They have also constrained the two 
parameters u>o and fi^o of the model using the SNLS 
(Supernva Legacy Survey) [2fjj and BAO datafUJ. At 
present, we have the Union08 compilation of the Snla 
data which contains around 307 data points [2^]. This is 
world's published first heterogeneous SN data set con- 
taining large sample of data from SNLS, Essence survey, 
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FIG. 5: Behavior of A as function of scale factor for different 
potentials. We have chosen the initial value of Xi = 0.5. The 
Dotted, dashed and dotdashed curves correspond to V(<fi) = 
(j) -1 , 4>~ 2 , (fi~ 3 respectively. 



high redshift supernova data from Hubble Space tele- 
scope as well as several small data sets. We use this 
data set together with the BAO data from SDSS (Sloan 
Digital Sky Survey) [21|. The la and 2a contour inter- 
vals for our model have been shown in figure 4. From 
the figure, it is clear that one can not distinguish cosmo- 
logical constant with a thawing dark energy models with 
present data although the phantom dark energy models 
are preferred. 



IV. CONCLUSIONS 

In this paper we have examined the DBI system with 
a phenomcnologically motivated class of run away po- 
tentials. In general, the the tachyon dynamics crucially 
depends upon the asymptotic behavior of the potential 
V{<f>) at large values of 4>. The inverse square potential 
gives rise to constant equation of state which is deter- 
mined by the slope of the potential, w = — 1 + A 2 /2. We 
analysed the class of tachyon potentials with dark energy 
and dark matter as late time attractors. Models in which 
V{4>) decrease faster than cj)~ 2 can give rise to transient 
dark energy near the top of the potential and then mimic 
dark matter as late time attractor. Since p^ for tachyon 
field scales slower than matter, its energy density for a 
viable cosmic evolution should be fixed around p\ at ear- 
lier epochs allowing the field to freeze due to large Hubble 
damping. Thus all the three classes of tachyon models 
belong to thawing type. The data available at present al- 
lows to carry out investigations around the present epoch 
with 70 << 1. As soon as p^ becomes, comparable to 



FIG. 6: The graph shows w-w' phase space occupied by the 
fields . The upper bound and the lower bound correspond to 
3(1 + w) and —0.8(1 + w) respectively. The curves are for 
the potentials V{4>) = (dotdashed curve), V(<p) = 4>~ 2 
(Short dashed curve), V(<j>) = (j)^ 1 (Dotted curve), V((f>) = 
(/r - 2 (Thick curve), V(<f>) = 0"° ^(Long Dashed curve).The 
two thin lines are represent the upper and lower bound for the 
thawing models. The corresponding bounds are also specified. 




FIG. 7: Constraints in wo — Jl^o parameter space using 
Union08 compilation of SN data and BAO data. Black lines 
are for SN data only while Red lines are for SN+BAO data. 
Solid lines are for la contour intervals while dashed lins are 
for 2cr contour intervals. 



matter density field begins evolve. The equation of state 
improves slightly starting from w(<f>) = — 1. Hence, the 
slope of the potential does not change appreciably which 
we confirmed numerically. In the limit of small adiabatic 
index of <p assuming A to be constant, we have shown that 
the resulting evolution equations are same as in case of 
quintessence which can be solved analytically. Our sim- 
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ulation shows that the approximation is very close to the 
numerical results for A < 1 around the present epoch. 
Deviations are possible in the far future. We therefore 
conclude that tachyon dynamics is difficult to distinguish 
from quintessence at least in the near future. We also 
extended our analysis to the case of phantom tachyon. 
Again in the region of interest, we find that phantom 
tachyon model is difficult to distinguish from the ordi- 
nary phantom field. We also constrained the parameters 
Wo and O^q for our model using the latest supernovae 
data along with baryon acoustic oscillation BAO data. 
Our analysis shows some preference for phantom energy. 

The fact that all the scalar field dark energy models 
have a unique equation of state as long as they are in 
the slow-roll regime, makes a strong case for the w(a) 
given by equations (17) and (19). It does not matter 



whether the scalar field has a canonical or non-canonical 
kinetic term. It is also the same for non-phantom or 
phantom scalar fields. We hope that this equation of 
state behaviour for the dark energy will be considered 
seriously while fitting with the observational data coming 
from future experiments. 
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